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Abstract 

We present the first measurement of the ratio of branching fractions R = 

5 



B(t Wb)/B(t -► Wg) from pp collisions at ^/s = 1.8 TeV. The data set 
corresponds to 109 pb _1 of data recorded by the Collider Detector at Fermi- 
lab during the 1992-95 Tevatron run. We measure R = 0.94+^24 (stat+syst) 
or R > 0.61 (0.56) at 90 (95)% CL, in agreement with the standard model pre- 
dictions. This measurement yields a limit of the Cabibbo-Kobayashi-Maskawa 
quark mixing matrix element \Vtb\ under the assumption of three generation 
unitarity. 

12.15.Ff, 12.15.Hh, 14.65.Ha 
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The Cabibbo-Kobayashi-Maskawa matrix |l] is a fundamental component of the stan- 
dard model of electroweak interactions. However, the matrix elements must be determined 
experimentally since the model does not constrain their values. Some of the matrix ele- 
ments have been determined by studying the weak decay of quarks or by deep inelastic 
neutrino scattering experiments. Until now no direct information has been available for the 
elements of the top sector. The matrix elements \Vtd\ and |Vt s | have been indirectly esti- 
mated and \ Vtb\ has been deduced by a global fit, with the additional assumptions of having 
only three generations and unitarity. With this procedure the indirect allowed range for 
|Vtb| is 0.9989 -T- 0.9993 (at 90% CL) 0. However, without the assumption of three quark 
generations, the constraint is far less stringent: \V t b\ — 0. -r- 0.9993 (at 90% CL) |2j. The 
large value of \V t b\ in the standard model implies that R, the ratio of branching fractions 
B(t —> Wb)/B(t — > Wq) (where q is a d, s or b quark), is close to unity and that the branch- 
ing ratio for the decay of a top quark to Wb is nearly 100%. This prediction has been used 
in the discovery of the top quark but has not been, until now, confirmed experimentally. 

In this letter we present the first direct measurement of R. The result provides additional 
support for the top quark discovery and the first direct constraint on the CKM element \V t b\ 
under the assumption of unitarity. The analysis is performed using 109 pb _1 of proton- 
antiproton collisions data recorded at a center of mass energy of 1.8 TeV by the Collider 
Detector at Fermilab (CDF) during the 1992-1995 run of the Tevatron Collider at Fermilab. 
The CDF detector is described in detail elsewhere ||; here we briefly describe only the 
components which play a major role in this analysis. 

The CDF tracking system consists of three different detectors embedded in a 1.4 T 
solenoidal magnetic field. A particle emerging from the interaction region passes through 
a four-layer silicon vertex detector (SVX) ||,|5| located just outside the beam pipe, a set 
of vertex time projection chambers (VTX) and a drift chamber (CTC) with 84 measuring 
planes. The CTC performs the pattern recognition and a three-dimensional reconstruction 
of charged particles. The VTX measures the position of the primary interaction vertex 
along the beam axis. Finally, the SVX, with its r — ip readout in the plane perpendicular 
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to the colliding beams, is designed to determine precisely the impact parameter of the 
tracks. A momentum resolution of APt/P? — 0.0009 (GeV/c) -1 || and an asymptotic 
impact parameter resolution of ~ 13 /xm is obtained for tracks with high Pt detected by the 
SVX and the CTC. Outside the tracking volume, electromagnetic and hadronic calorimeters 
measure the energy of particles in the region < 4.2 J7|. Electron identification is obtained 
by combining calorimetric and tracking information. Muon identification is performed by 
matching tracks reconstructed in the CTC with segments measured by a system of drift 
chambers (muon chambers) located outside the calorimeters and covering the region |?7| < 1. 

The top quark has been observed |]8|-[T0l only when produced in pairs. Assuming that the 
top quark decays to a real W boson, it is customary to classify tt final states according to 
the decay modes of the two W bosons. We use two tt candidate data sets: the lepton+jets 
(l+jets) and the dilepton samples. The l+jets sample, in which one W decays to an electron 
or a muon and its corresponding neutrino and the other W decays into two jets, has a final 



state characterized by a high-P^ lepton, missing transverse energy Ex |jl 1|| 7 and four jets. 
The dilepton sample, in which both W bosons decay into electrons or muons and neutrinos, 
is characterized by a final state with missing transverse energy, two high-P^ leptons and two 
jets. The selection criteria required to identify the W bosons and to enhance the top quark 
content in these data sets have been described in detail in previous CDF publications |||| • 
To isolate the I + jets sample we require the presence of one central (\rj\ < 1) lepton (e or 
H) with P T > 20 GeV/c, E T > 20 GeV, three jets with E T > 15 @ GeV within \rj\ < 2 
and a fourth jet with > 8 GeV within < 2.4. In the dilepton sample we require two 
leptons (e or //) with P T > 20 GeV/c, E T > 25 GeV and two jets with E T > 10 GeV in the 
region \rj\ < 2.0. Candidate Z events are removed by rejecting events containing same-flavor 
lepton pairs with opposite charge whose invariant mass lies between 75 and 105 GeV/c 2 . By 
construction the two data sets have no overlap. After applying all the selection criteria we 
find 163 events in the I + jets and 9 events in the dilepton sample. 

The presence of the t Wb decay is deduced by identifying jets associated with b hadron 
decays using two distinct algorithms: the SVX tagger and the SLT tagger. The SVX tagging 
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algorithm || relies on the long lifetime of b hadrons. It searches for b hadron decay vertices 
which are significantly displaced from the primary vertex and have three or more associated 
tracks. If this search fails, tighter quality cuts are applied to the tracks within the jet and 
vertices with two tracks are also accepted. In both cases, the transverse displacement of the 
decay vertex from the primary vertex, divided by its uncertainty, is required to be larger 
than 3. The SVX algorithm is characterized by an efficiency (e s ) to tag a single b jet in 
a tt event of (37.0 ±3.7)% and by a fake tagging rate of about 0.5%. The SLT tagging is 
performed by looking for 1ow-Pt (relatively soft compared to the primary lepton) muons 
and electrons from semileptonic b hadron decays. The algorithm looks for low transverse 
momentum electron and muon candidates by matching CTC tracks with > 2 GeV/c 
with calorimeter clusters and track segments in the muon chambers. Moreover, to classify 
the event as a tt candidate, the soft lepton is required to be within AR < 0.4 |13| from one 
of the four highest-Er jets in the event. The SLT algorithm has an efficiency per jet (ei) of 
(10.2 ± 1.0)% and a fake tagging rate of about 2%. Background due to fake tags is measured 



for both algorithms using samples of QCD jet data sets \TA . 

The unknown ratio of branching fractions, R, is measured by comparing the observed 
number of tags in the data with expectations based on selection criteria acceptances, tagging 
efficiencies and background estimates. In the dilepton sample only SVX tagging is used and 
the sample is divided into three non-overlapping bins: events with no 6-tags (bin 0), one 
and only one 6-tag (bin 1) and two 6-tags (bin 2). The use of SLT tagging in the dilepton 
data set does not provide any additional statistical gain. In the / + jets sample, we use 



both the SVX and SLT algorithms. Monte Carlo studies |15|] indicate that a superior use of 
the tagging information is obtained by dividing events into the same three bins used for the 
dilepton sample and then by subdividing the bin with no SVX tags into two bins according 
to the SLT tagging status. The first bin (bin 00) is populated by events which are tagged 
by neither the SVX nor the SLT algorithm and the second one (bin 01) contains events with 
one or more SLT tags and no SVX tags. 

The number of observed events in each bin is reported in Table |. The expected number 



of events, iVj, in each of the bins of the I + jets sample can be expressed as a function of 
the acceptances, tagging efficiencies and the estimated background, by the following set of 
equations: 

iV 00 = n + (1 - e,)(l - e.)ni + (1 - e{) 2 {l - e s ) 2 n 2 + F 0Q (la) 
iV i = ei{\ - e 8 )m + ei{2 - e,)(l - e s fn 2 + F m (lb) 
Nx =e s n 1 + 2e s {l-e s )n 2 + F 1 (lc) 
iV 2 = ^ n2 + F 2 ( ld ) 

with rii (i = 0, 1, 2), the number of events with % 6-jets in the SVX acceptance, given by: 

n = N top [a + (1 - R)a! + (1 - R) 2 a 2 ] (2a) 
^ = N top [R ai + 2/2(1 - i?)a 2 ] (2b) 
n 2 = N top R 2 a 2 (2c) 

where iV to p is the total number of tt events in the sample, Fi is the background in the i-th 
bin and a, is the fraction of events containing % 6-jets (i = 0, 1,2) in the acceptance. This 
definition of acceptance, which reflects the way the c^'s are related to R in Eq. (fZa|)-(|2"c|), 
has been chosen in order to be able to use the standard CDF top Monte Carlo (see below) 
which assumes R = 1. For the dilepton sample, Eq. ([la]) and (|Tb|) are merged into one 
because SLT tagging is not used. 

The unknown ratio R is obtained by minimizing the negative logarithm of a likelihood 
function. Since the I + jets and dilepton samples are independent, the global likelihood can 
be written as: 

C> = £i+jets£ dilepton (3) 

where each of the individual likelihoods is of the form: 

C = ]Jl'i\r.\ i )l[(;(.r J :.r r a J ) 
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In this expression, P(Ni\ N) is the Poisson probability for observing iVj events in each 
bin (the index i runs from 1 to 4 for the I + jets sample and from 1 to 3 for the dilepton one) 
with an expected mean Ni (see Table |). The functions G(xj] Xj,aj) are Gaussians in Xj, 
with mean Xj and variance crj, and incorporate the uncertainties in the tagging efficiencies, 
backgrounds and acceptances into the likelihood functions. 

The acceptances and efficiencies are obtained using a ti Monte Carlo (M top = 175 
GeV/c 2 ) data set generated using PYTHIA [li|, combined with a detailed simulation of 



the detector response. The total number of ti pairs (N top ) in the two data samples is left as 
a free parameter. The acceptances in each bin are normalized with respect to the bin with 
no 6-jets. As a consequence, the trigger and lepton identification efficiencies cancel out in 
the ratio. We obtain r x = 11.8 ± 1.2 (14.5 ± 1.4) and r 2 = 38.7 ± 3.9 (58.5 ± 5.8), where 
ti = cii/cio, for the I + jets (dilepton) sample. The uncertainties in these ratios include 
contributions from the jet energy scale and from the Monte Carlo modeling of initial and 
final state radiation. 

The background in the untagged sample is mainly due to the associated production 
of W bosons with light quark jets. The backgrounds to the SLT and SVX tagged events 
(background in bin 01 and 1, respectively), are mainly due to the associated production of 
W bosons and heavy quarks (Wbb, Wcc, Wc) and to mistags due to mismeasured tracks. 
Smaller contributions come from bb, diboson production (WW, ZZ and WZ), Z — > rr de- 
cays, Drell-Yan lepton pair production and single top quark production. These backgrounds 



are calculated using a combination of data and Monte Carlo information [p^ , [15|| . The initial 
values of the SVX and SLT backgrounds are a function of the ti content of the / + jets 
sample itself, and therefore need to be appropriately corrected ||. An iterative process is 
used to account for this effect and has been implemented in the likelihood minimization pro- 
cedure used to estimate R. Using this procedure, as output of the likelihood minimization, 
we estimate F% = 3.3^2 and F i = 7.2 ± 1.6 events for the SVX and SLT backgrounds, 
respectively. In the same way, the background to double SVX tagged events (bin 2) is 
estimated to be small and amounts to F2 = 0.2 ± 0.1 events. The background in bin 00, 
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Foo, is obtained as follows. Defining N tot to be the total number of events in the I + jets 
data set, N SV x the total number of SVX tagged events in this sample, F S vx the estimated 
background and e$vx the SVX event tagging efficiency, the total number of top events is: 
N t t = (Nsvx — Fsvx) / (tsvxR)- The number of background events before tagging is given 
by: F = N tot — N t i and therefore the background in bin 00 is F 00 = F — (F 01 + F\ + F 2 ). As 
before, the estimate is performed iteratively during the likelihood minimization. 



The initial background to the dilepton sample has been estimated to be 2.4 ± 0.5 |L7 
In this sample, we estimate a background of 0.10 ±0.04 events to SVX single tagged events. 
The double SVX tagged background is negligible (F 2 = in Eq. ([id])). In this case, the 
number of background events is not a function of the tt content of the initial sample and 
no special correction need to be applied. As in the / + jets case, the background in bin 
is obtained by a subtraction of the tagged background from the total background and 
amounts to 2.3 ± 0.5 events. The resulting number of background events after the likelihood 
minimization procedure is shown in Table |TJ for both data sets. 

The likelihood minimization yields R = 0.94+Q24. The uncertainties includes both statis- 
tical and systematic effects with the former being the dominant contribution. The negative 
log-likelihood as a function of R is shown in Fig. The lower limit on R is obtained by a 
numerical integration of the likelihood function and we obtain: R > 0.61 (0.56) at 90 (95)% 
CL. 

The CKM element \V t b\ is directly related to R, although in a mo del- dependent way. We 
assume that the top quark decays to non-W final states are negligible |T8|JT9|| . Under this 
assumption R is related to \V tb \ by: 





V* 


to 


\V U 


2 + 


V td 


2 + 


v tb 


2 



If we assume three generation unitarity, the denominator is equal to unity and therefore 
R = \V tb \ 2 . As a consequence, we obtain: \V tb \ = O.Q7t°it or \V tb \ > 0.78 (0.75) at 90 (95)% 
CL. 

The result, although limited by statistics, represents the first direct measurement of 
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R. The large value of R that we measure is consistent with standard model expectations 
and supports the assumption that top quarks decay predominantly to b quarks. Under the 
assumption of three generation unitarity, our calculated value of | V tb | = 0.97t^ (\V tb \ > 0.78 
at 90% CL) is consistent with indirect limits obtained from global fits. 
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TABLES 



sample 


00 


01 




1 2 


I + jets 
dilepton 


126 
6 


14 

n/a 




18 5 
3 


TABLE I. Number of events for the two data samples in each of the bins defined according to 
the number of SVX and SLT tags. The case of SLT tags (bin 01= "one or more SLT tag, no SVX 
tag") does not apply ("n/a") to the dilepton data set (see text). 


sample 


00 


01 




1 2 


I + jets 
dilepton 


108 ± 10 
2.3 ±0.5 


7.2 ± 1.6 

n/a 




3.3tH 0.2 ±0.1 
0.10 ±0.04 n/a 



TABLE II. Estimated number of background events bin by bin for the two data samples. The 
case of SLT tags (bin 01) applies only to the I + jets data set and the double tagged dilepton 
background is neglected ("n/a"). 
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